Lateral transport of thermal capillary waves 
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We demonstrate that collective motion of interfacial fluctuations can occur at the interface between 
two coexisting thermodynamic phases. Based on computer simulation results for driven diffusive 
Ising and Blume-Capel models, we conjecture that the thermal capillary waves at a planar interface 
travel along the interface if the lateral order parameter current j{y) is an odd function of the distance 
y from the interface and hence possesses opposite directions in the two phases. Such motion does not 
occur if is an even function of y. A discrete Gaussian interface model with effective dynamics 
exhibits similiar transport phenomena but with a simpler dispersion relation. These findings open 
up avenues for controlled interfacial transport on the nanoscale. 

PACS numbers: 05.40.-a, 05.50.+q, 68.05.Cf, 68.35.Rh 



Understanding the motion of interfaces is important 
in areas such as multiphase flow, dendritic and crystal 
growth, microchip fabrication, combustion, blood flow, 
and cell dynamics. Pronounced changes to the micro- 
scopic interfacial structure can occur as a result of the 
motion, e.g., the development of kinetically enhanced 
self-afSne roughness [ij of the driven interface between 
a growing material and its environment. In striking con- 
trast, driving the interface along the plane of its average 
position can lead to interfacial smoothening, as has been 
observed experimentally for a colloidal gas-liquid inter- 
face under shear flow and in computer simulations 
of interfaces in driven lattice gas models d, |3| • Coher- 
ent lateral transport of large-scale interfacial structures 
is well-known in macroscopic systems that are far from 
equilibrium, such as migrating sand dunes and ripples, 
or ocean waves. Finding thermal analogues can be of rel- 
evance in micro- and nanofluidic devices, e.g. due to the 
possibility of transport of nanoparticles at a liquid-liquid 
interface by controlled motion of thermal capillary waves. 

Consider two distinct equilibrium phases at thermo- 
dynamic phase coexistence, organised into two different 
regions of space that are separated by a planar interface. 
The local order parameter $ varies along the axis y per- 
pendicular to the interface plane and for large enough 
systems reaches its respective bulk value far from the in- 
terface. The (scalar) order parameter profile $(?/) can 
represent the density near a gas-liquid interface, the rel- 
ative concentration at the interface separating two coex- 
isting liquid phases, or the magnetization near an Ising 
domain wall. On a coarse-grained scale, the interface can 
be characterized by its local departure (height) h{x,t) 
from a reference plane y — h{x, t) = 0, where x indicates 
the coordinate(s) parallel to the interface and t is time. 



At temperatures T above the roughening transition, the 
interface exhibits large spatial fluctuations. Phenomeno- 
logical capillary wave theory Q, as well as rigorous re- 
sults 0, indicate that in the absence of external fields 
that couple to the order parameter, the length scale of 
such fluctuations diverges with system size, i.e., the in- 
terface thickness becomes infinite in the thermodynamic 
limit in spatial dimensions d < 3. 

In this Letter we consider external driving that cre- 
ates a steady state with non-vanishing current of the or- 
der parameter, j{y), parallel to a planar interface. We 
investigate the effects on the dynamics of the capillary 
waves via the two-point correlation function C{x, t) = 
(/i(0, Q)h{x, t)), where the angles denote an average in the 
steady state. Based on computer simulations results for 
C{x,t) for various simple microscopic models, we con- 
jecture that the lateral flux of the order parameter at 
a planar interface induces lateral motion of the thermal 
capillary waves, provided j{y) is an odd function of the 
distance y from the interface. The spatial symmetry at 
equal times, C{x, 0) = C{—x, 0), is then broken for times 
t > 0, such that C{x,t) ^ C{—x,t). For general forms 
of j{y) it is the odd component, [j{y) — j{—y)]/2, that 
induces the motion. No motion occurs when jly) is even. 
We shed light on the transport mechanism by construct- 
ing a corresponding (Gaussian) effective interface model 
that displays similar transport phenomena. 

As a microscopic approach we use kinetic lattice mod- 
els with dynamics that conserve the order parame- 
ter locally. The Ising and the Blume-Capel [7|] mod- 
els of binary mixtures possess the Hamiltonian Ti = 
—Jj2{i I) "■jO'ij where (i, I) denotes nearest-neighbor pairs 
and the spin-spin coupling constant J > 0. Contribu- 
tions to H from static external fields are omitted, as we 
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will work in an ensemble of fixed numbers of spins cr^ 
of each type. For the Ising model ai = ±1, while for 
the Blume-Capel model cji — — 1,0,+1, which in lat- 
tice fluid language corresponds to occupancy of site i 
by particles of species "—1", vacancy, or species "+1" 
of the mixture, respectively. We initially consider two- 
dimensional square lattices of dimensions x Ly with 
coordinates x (horizontal) and y (vertical). The inter- 
face with mean orientation in the x-direction is estab- 
lished and localized by assuming boundary conditions 
CTj = -1-1 at the upper {y = {Ly + l)/2) and ai = —1 
at the lower {y = —{Ly + l)/2) edges of the lattice. Peri- 
odic boundary conditions are applied in the x direction. 
As a consequence, $(?/) = Uij/L^, i = {x,y) crosses 
over from <i> < for y < to <& > for y > 0. We 
use Kawasaki spin exchange dynamics Q , where the ele- 
mentary move consists of swapping the values of the spin 
variables Ui and cr; of two nearest-neighbor sites i and I. 
The system is driven by a force field F{y) that acts in 
the a;-direction (parallel to the interface) and varies with 
distance y from the interface. The acceptance rates for 
the trial moves are assumed to be of modified Metropolis 
type, min{l,exp(-(A7? + I^W) / {ksT))}, where Ai? is 
the change in internal energy, /S.W is the work due to the 
driving field, and ks is the Boltzmann constant. Hence 
the work performed by the field is dissipated into a heat 
bath, which is kept at T = const. We have carried out ex- 
tensive Monte Carlo (MC) simulations using multi-spin 
coding techniques extended for Kawasaki dynamics 
to include drive. Results are presented for L^ = 200 
and Ly ~ 20 at fixed total magnetization Ui — Q for 
T/Tc = 0.75 (where = 2.2692 J//cb) and run lengths of 
the order of ^Vmc = lO'^ MC steps {L^ x Ly trial moves). 

We first discuss the Ising lattice gas where the drive 
creates a work term for an exchange in the x-direction, 
AW = -JF{y){(Ti - ai)/2, where i = {x,y) and I = 
(a;-|- 1, y); exchanges in the y-direction occur with normal 
equilibrium rates, AW — 0. In the case of odd symmetry 
upon spatial reflection, F{y) = — F(— y), the field acts in 
opposing directions in both halves of the system, e.g. 
with linear variation across the slit, F{y) = jy, where 
7 is the (scaled) field difference between rows. Simula- 
tion results indicate that the order parameter current 
profile j(y) = j+(y) - j_(y), where ja{y) is the par- 
ticle current profile of species a ~ ±1, possesses the 
same direction and the same symmetry as the driving 
field. For odd driving fields, we find j{y) — —j{—y)- We 
define the instantaneous interface position via a coarse- 
graining method based on the (scaled) column magneti- 
zation, h{x,t) — —{2'mh)~^^yCFi, where site i — {x,y) 
possesses the value Ui at time i, and mb is the sponta- 
neous equilibrium magnetization in bulk. 

In Fig. [1^ we plot C{x, t) as a function of x for several 
fixed values of time difference t for the case of linear drive, 
7 = 1. At t = 0, a cusp is apparent at a; = and C{x, 0) 
exhibits long-ranged decay with distance x, characteriz- 
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FIG. 1: (color online) a) Height-height correlation function 
C{x,t) of the Ising model as a function of the distance x 
at constant time t (as indicated, in units of MC steps) for 
temperature T/Tc = 0.75, and system size Lx = 200, and 
Ly — 20. The system is subject to linearly varying drive with 
strength 7 = 1. Inset: C{x,t) for the case of V-shaped drive 
for the same parameters, b) Same as a), but for the discrete 
Gaussian interface model subject to linear /i-dependent drive. 
Inset: Illustration of the motion of the fluctuating interface. 



ing the (equal-time) spatial correlations of the interfacial 
fluctuations. Upon increasing time, t > 0, the position of 
the peak moves towards negative values of x, its height 
decays and its width increases. This behaviour clearly in- 
dicates the existence of damped propagating modes that 
move in the negative x-direction. The position of the 
maximum of C{x,t) varies linearly with time. The in- 
ferred velocity is Wpeak = 0.009 (in units of lattice con- 
stant per MC step); for not too large values of 7 the ve- 
locity f peak grows linearly with increasing 7. For step-like 
drive of strength /, i.e., F{y) = /sgn(y), where sgn(-) is 
the sign function, we observe very similar behaviour of 
C{x,t) (data not shown), from which we conclude that 
the occurrence of the interfacial motion is not tied to the 
specific functional form of the (odd) driving field. Quali- 
tatively different behaviour occurs for even symmetry of 
the drive, F{y) = F{—y), such that the drive acts in the 
same direction throughout the system. We find, for the 
cases considered, that the order parameter current pro- 
file also attains even symmetry, j(y) = j{—y). For a V- 
shaped spatial dependence, -F(y) = 7|y|, simulation data. 
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FIG. 2: (color online) Dispersion relation for the lat- 

eral propagation of capillary waves in the two-dimensional 
Ising model driven by a linear field for 7 = 1 (a), and in the 
one- dimensional kinetic discrete Gaussian model with height- 
dependent driving (b). The spectral variable is g = 2-nk/Lx, 
where k is the wave number. Parameters are the same as in 
Fig. [1] Data for different system sizes = 100, 200, 400 col- 
lapse onto each other. Also shown is the analytical formula 
(line) described in the text. 



FIG. 3: (color online) a) Density profiles Pcr{y) of species — 1, 
-1-1, and vacancies (0) for the Blume-Capel model under ho- 
mogeneous drive of strength / = 0.25 and T = 0.75J/fcs 
as a function of the scaled distance from the mean inter- 
face position y/Ly. Also shown is the order parameter profile 
^(y) ~ P+{y) ~ P-iy)- inset: C(x,t) as a function of x for 
various values of t (as indicated), b) Particle current profiles 
jaiy) for species a = —1,-1-1,0 and order parameter current 
profile, j{y) = j+{y) - j- (y), as a function of y. 



shown in the inset of Fig. [T^, show that with increasing 
time the peak of C{x,t) decreases in magnitude but re- 
mains stationary at x = 0, which indicates the absence of 
propagating modes. This observation holds also for the 
case of uniform drive, F{y) = / = const. The intercept of 
the equal-time correlation function, C(0,0) = (/i(0,0)^, 
provides a measure of the (squared) interfacial width [J| ; 
we find for all cases considered that C(0,0) decreases 
under drive. Comparison of the results for linear and V- 
shaped drive indicates that the suppression of roughness 
is less strong in cases where interfacial motion occurs. 
The interfacial transport is intimately related to a broken 
symmetry under space reflection, y — > —y and x — > —x, 
and species inversion, (Ti ~ai, that occurs for cases 
of odd driving. As a consequence, on average —1 spins 
in the region with y > move in the same (negative) x 
direction as -f 1 spins do in the region with y < 0. Hence 
"intruders" of the other phase move in the same direc- 
tion throughout the system. This direction is opposite 
to that of the "velocity profile" of the order parameter, 
Vit>{y) = j{y)l^{y)^ which is O < v^ijj) < 0.022 for all y. 
See the inset of Fig. for an illustration of the transport 



phenomenon. 

The discrete Gaussian model [l^ provides a reduced 
description with interfacial degrees of freedom only. Nev- 
ertheless these are known to exhibit long-wavelength 
equilibrium fluctuations that are characteristic of real 
fluid interfaces. In a one-dimensional system the inter- 
face is represented by integer height variables h{x) with 
X — \,...^Lx and periodic boundary conditions. The 
Hamiltonian is T^dg — {J /'^)Yli^=i{^{^ ^ ^) " ■ 
In our (conserved) dynamics a site x and one of its neigh- 
bours x' ~ X ±1 are chosen at random, and the heights 
are changed as h{x) h{x) + 1 and h{x') — > h{x') — 1, 
with probability given by the modified Metropolis rate. 
This can be viewed as an (effective) biased diffusion step 
of a particle along the interface from x' to x. We have 
performed MC simulations for the same parameters as for 
the Ising model, and for two cases of driving. In the first 
case, the work is AW = -J^{x' - x){h{x) + h{x'))/2, 
corresponding to linear drive in the Ising lattice gas. 
As shown in Fig. [Hd, the temporal height-height cor- 
relation functions exhibit characteristics very similar to 
those found in the microscopic (Ising) model with cap- 
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illary wave motion. The velocity of the position of the 
maximum is Wpcak — 0.08 (in units of lattice constants 
per MC step, where one MC step consists of exchange 
moves) , significantly larger than that of the Ising model. 
Biasing the moves to the left (say) irrespective of the 
height variable, AW = ~'J{x' — x)f, corresponds to uni- 
form drive in the Ising model. As in the Ising model, no 
transport is observed. In cases where transport occurs, 
the combined symmetry h —h, x —x, and exchange 
of X and x' , is broken. 

The short-time dynamics of the capillary waves can 
be characterized by a dispersion relation of frequency uj 
as a function wave number q = 2Trk/Lx, obtained from 
the (average) phase shift of mode k = 1,2, . . . , in 
unit time as w{q) = a,Tg{h* {q,t)h{q,t -\- dt))/dt, where 
h{q,t) is the spatial Fourier transform of h(x,t), and dt 
is a small time interval (L^j/lO spin exchanges). Sur- 
prisingly, for small wavenumbers q the driven Ising lat- 
tice possesses a nonlinear dispersion relation, see Fig. [2^ 
where results are shown for the case of linear drive with 
7 = 1. In contrast, uj{q) for the discrete Gaussian model 
is linear for small q, see Fig. ^p. Modelling the dynam- 
ics by a simple linear transport operator dt — vdx, with 
the continuous partial time derivative dt, discrete spa- 
tial derivative dxh{x,t) = [h{x + l,t) — h{x — l,t)]/2, 
and plane wave modes exp(i(wt + qx)) yields uj{q) = 
vsm{q). This describes the simulation data very well, 
see Fig. [2j). For small q, the behaviour uj{q) — vq, 
with V = 0.0760(1), is in good agreement with the 
value of Wpoak obtained from analysis of C{x,t). The 
much richer behaviour of the Ising model can be fit- 
ted to u){q) = {v -\- 2u) sui{q) — usin(2g) -|- ssin^(q) with 
small-(? expansion ijj{q) = vq + sq^ + {u — v/6)q^, and 
V = 0.0172(7), ti = 0.0402(2), s = 0.0263(13). A linear 
transport operator dt — vdx +udx (5-point stencil for 9^) 
yields the first and second terms of the fit function, but 
the third term cannot be obtained from a linear transport 
equation with real coefficients [ll[. We expect a quanti- 
tative description of the full dynamics, as measured by 
C{x,t), to require a non-linear description. 

In the kinetic Blumc-Capel model the presence of va- 
cancies (sites i with ai = 0) creates richer microscopic 
dynamics. In contrast to the Ising model, where driv- 
ing -|-1 particles is intimately related to counter-driving 
— 1 particles in the opposite direction, the Blume-Capel 
model offers the possibility to co-drive the ±1 species 
in the same direction. The corresponding driving fields 
F„{y) acting on species ct = ±1 then obey F+{y) — 
F^(y) = F{y). The work term for an exchange of spins 
i and I with coordinates x and x' = x + 1, respectively, 
is AW = —JF{y){a^ — af). Simulation results for va- 
cancy concentration of 20%, where we find the confined 
interface to be stable, indicate that the symmetries of the 
driving field F{y) and of j(y) are no longer the same, e.g., 
that an even field can give rise to an odd order param- 
eter current profile, see Fig. [3] for results of the density 



profile, Pa{y), and current profile, ja{y), of each species 
a = —1, 0, -1-1 for the case of uniform drive, F{y) = const. 
Indeed motion of capillary waves is observed, see the in- 
set of Fig. [3K for results for C{x, t). We do not find wave 
motion when counter-driving, i.e. for F^[y) = —F-(y). 
All these findings are fully consistent with the proposed 
scenario that the interfacial motion is caused by the sym- 
metry of j{y) and not by that of the driving fields. 

For general cases where the driving field and hence the 
order parameter current profile do not possess unique 
but rather mixed spatial symmetry, we find, for the Ising 
model, that it is the presence of an odd component that 
leads to interfacial motion. We have also checked that the 
transport phenomenon is not specific to two-dimensional 
systems, but also occurs in the three dimensional Ising 
model, where simulation results indicate that the direc- 
tion of motion is parallel to that of the drive and that the 
second coordinate parallel to the interface plays a mere 
spectator role. 

In view of possible microfluidic applications, it would 
be very interesting to test the validity of the proposed 
scenario in experiments, e.g., using colloidal dispersions 
2], and in molecular dynamics simulations [l^ . 
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